Nonclassicality and decoherence of photon-subtracted squeezed states 
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We discuss nonclassical properties of single-photon subtracted squeezed vacuum states in terms 
of the sub-Poissonian statistics and the negativity of the Wigner function. We derive a compact 
expression for the Wigner function from which we find the region of phase space where Wigner 
function is negative. We find an upper bound on the squeezing parameter for the state to exhibit 
sub-Poissonian statistics. We then study the effect of decoherence on the single-photon subtracted 
squeezed states. We present results for two different models of decoherence, viz. amplitude decay 
model and the phase diffusion model. In each case we give analytical results for the time evolution 
of the state. We discuss the loss of nonclassicality as a result of decoherence. We show through the 
study of their phase-space properties how these states decay to vacuum due to the decay of photons. 
We show that phase damping leads to very slow decoherence than the photon-number decay. 

PACS numbers: 03.65.Yz,42.50.Dv 



I. INTRODUCTION 



Quantum states can be well described in terms of 
Wigner functions. The states with Gaussian Wigner 
function have been of particular interests in context of 
quantum information processing. Entanglement in such 
states in terms of quadratures are also well studied. On 
the other hand, the quantum states with non-Gaussian 
Wigner function are also quite important. For exam- 
ple, a single-photon state, which finds many applications 
in quantum information processing, shows non-Gaussian 
behavior in phase-space. In a recent experiment, a non- 
Gaussian state has been produced by homodyne detec- 
tion technique from a single-mode squeezed state of light 
[l], 0. For certain non-Gaussian states, Wigner func- 
tions can take negative values. Such negativity refers to 
nonclassicality of these states. These states are useful 
in entanglement distillation 0, 0] , loophole- free tests of 
Bell's inequality @, and quantum computing [6[. A spe- 
cific class of such nonclassical states has been shown to be 
similar to the Schrodingcr kitten state, in the sense that 
their Wigner functions show negativity at the origin of 
phase space @, H| ■ It is well known that the Schrodinger 
cat states @, which are quantum superpositions of co- 
herent states, are non-classical in nature and are very 
important to study the interface of quantum and clas- 
sical worlds. Superposition of coherent states with low 
amplitudes creates Schrodingcr kitten states. Most of the 
experiments to prepare the Schrodinger cat states have 
been performed in cavities or bound systems. Thus they 
are not much useful in quantum information networks 
though they have the non-Gaussian nature which is re- 
quired in certain quantum communication protocols. In 
0, Q , it has been shown how to prepare an Schrodinger 
kitten state in an optical system, by subtracting a sin- 
gle photon from a squeezed vacuum state. This optical 
kitten state would overcome the limitations of bound sys- 
tems. Repeated photon-subtractions can lead to condi- 
tional generation of arbitrary single- mode state [l(| • We 



note that similar non-Gaussian states could be prepared 
by adding a single photon to a squeezed vacuum state 
(see [ill, G3 for details of photon-added coherent states, 
which are also non-Gaussian states). These states are 
equivalent to single-photon subtracted squeezed vacuum 
state and exhibit similar behavior in phase space. It is 
worth to mention that non-Gaussian two-mode entangled 
states can be prepared by subtracting a photon from a 
two-mode squeezed state [H, . 



In this paper, we focus our study on the nonclassical 
properties and decoherence of single-photon subtracted 
squeezed vacuum states which are optically produced 
single-mode non-Gaussian states. The structure of the 
paper is as follows. In Sec. II, we introduce the photon- 
subtracted squeezed states and discuss its nonclassical 
properties in terms of the sub-Poissonian statistics and 
the negativity of its Wigner function. We derive a com- 
pact expression of the Wigner function and find the re- 
gion in phase space where it becomes negative. We show 
that there is an upper bound of the squeezing param- 
eter for this state to exhibit sub-Poissonian statistics. 
In Sec. Ill, we study the effects of two different model 
of decoherence: photon-number decay and phase damp- 
ing. In both cases, we derive analytical expressions for 
the time-evolution of the state and its Wigner function. 
We discuss the loss of nonclassicality due to decoher- 
ence. We show through the study of evolution of the 
Wigner function how the state decays to vacuum as a 
result of photon-number decay. We further show that 
phase damping leads to much slower decoherence than 
the photon-number decay. 
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FIG. 1: (a) Gaussian Wigner function W sq of the squeezed vacuum state \ip) for 8 = and r = 0.31. (b) The Gaussian Wigner 
function of a vacuum state. 

FIG. 2: Plots of Wigner functions of single-photon subtracted squeezed states for (a) r = 0.31 and (b) r = 0.8 with 8 = 0. 



FIG. 3: Variation of C in phase space for r = 0.31 and 8 = 0. FIG. 6: Variation of P/(a 2 - 4|c| 2 ) 5/2 in phase space for r = 

0.31, 8 = 0, and Kt = 0.1. 

FIG. 4: Contour plot for C = constant in phase space for (a) 

r — o 3i anc j ^ r — o g FIG. 7: Plot of 16C' /(l — e~ 2Kt ) in phase space for (a) r = 0.31 

and (b) r = 0.8 at times (i) Kt = 0.05, (ii) Kt = 0.1, (iii) 
nt = 0.3, and (iv) nt = 0.5. 



II. PHOTON-SUBTRACTED SQUEEZED 
STATES 



An unnormalizcd single-mode squeezed vacuum state 
is given by 



\i;)=S(0\0) , S(C)=exp 



c 

r 



f2 



(1) 



where 5(C) is the squeezing operator, £ = re 1 ® is the 
complex squeezing parameter, and a is the annihilation 
operator. The Wigner function of this state is Gaussian 
and positive in phase space. When p (> 0) number of 
photons are subtracted from such states, the state can 
be written as 



|0> 



(2) 



where £ = tanh(£). For odd p — 2m + 1, the normalized 
form of this state can be written as 



(C/2) 



(2s + 2m + 2) 



N j^ (s + m + 1)! V2iTT 
while for even p = 2m, the state becomes 
1 ^ (^/2) fl + m (2s + 2m)! 



I2s + 1) 



2s) 



(3) 



(4) 



where N a and iV e are the normalization constants. In 
this paper, we focus on the case when a single photon is 
subtracted from the squeezed vacuum, i.e., for to = in 
©■ 



A. Negativity of the Wigner function 
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FIG. 5: Variation of the Q-parameter with time in presence of 
decoherence due to decay of photon for squeezing parameters 
r = 0.31 (red line) and r = 0.8 (green line). 



The Wigner function of the squeezed vacuum state \ip) 
is given by 

W sq (a ia *) = -exp(-2|a| 2 ) , (5) 

where a — acosh(r) — a*e sinh(r). The function is 
Gaussian in phase space, as shown in Fig. [TJ We now 
calculate the Wigner function of the state \if)) p - This 
state can be rewritten as 

\i/>) p = a p S(C)|0) = S(QSt(Qa?S(£)\0) . (6) 
Using the relation 

S*t (Oa p S(() = [cosh(r)a + e w sinh(r)a t ] p , (7) 
we get the following from © 

Wi=S(C)|l>, (8) 
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FIG. 8: Wigner function of single-photon subtracted squeezed states for 9 = and r = 0.31 at (a) nt = 0.05, (b) Kt = 0.1, (c) 
Kt = 0.3, and (d) Kt = 0.5. 



FIG. 9: Wigner function of single-photon squeezed states for 9 — and r = 0.8 at (a) Kt = 0.1, (b) Kt = 0.3, (c) Kt = 0.5, and 
(d) Kt = 0.7. 



for p = 1. For a density matrix p(a, cr), we can write the 
following: 

S(C)p(o, at)St(C) = p[5(C)a5t(C), 5(0^(0] . (9) 

We write the Wigner function of p(a,a^) as Wp(cv,a*). 
Using the identities 

S(QaS*(() = acosh(r) -aV e smh(r) , (10) 
£(00*5* (0 = a f cosh(r) - ae" 19 sinh(r) , (11) 

in ([5]) , we can thus write the Wigner function of the den- 
sity matrix p{a,a { ) — S{C)p{a,a [ )S^ {Q as 



W p (a,a*) = Wp(a,a*) 



(12) 



where we have used the linearity property of the Wigner 
function. For the state (JSJ), p = |1)(1| and its Wigner 
function is given by 



Wp(a,a*) = -(4|a| 2 - l)e 



-2|c 



(13) 



Thus, the Wigner function of the single-photon sub- 
tracted squeezed vacuum state becomes 



W p (a,a*) = -(A\a\ 2 -l)e 



-2|«| 2 



(14) 



where we have used (|T2|) and (|13j) . Clearly, the Wigner 
function (|14p is non-Gaussian in phase-space. We show 
the plot of this Wigner function in the phase-space in 
Figs. [2] for different squeezing parameters. As an evi- 
dence of non-classicality of the state, squeezing in one of 
the quadratures is clear in the plots. Also there is some 
negative region of the Wigner function in the phase-space 
which is another evidence of the non-classicality of the 
state. The function becomes negative in phase space, 
when 

' (15) 



|a| 2 < 



U,|2 



We show in Fig. [3] the variation of C 
phase space. The negative region of C corresponds to 
the negativity of the Wigner function. Note that C = 
constant corresponds to ellipse in phase space, as shown 
in Fig. |U 

Note that the photon-subtracted squeezed states are 
similar to Schrodinger kitten states @, HI because Wigner 
functions exhibit the same characteristics in phase-space 
especially for the large values of the squeezing parame- 
ters. Moreover, in both cases, Wigner function becomes 
negative in the center of phase space. 



B. Sub-Poissonian nature of the photon-subtracted 

state 

The nonclassicality of the state can also be ana- 
lyzed by studying its sub-Poissonian character in terms 
of the Mandel's Q-parameter [15[ which is defined by 



Q 



(oto) 



(16) 



The negativity of the Q-parameter refers to sub- 
Poissonian statistics of the state. However in 16], it has 
been shown that a state can be nonclassical even if Q is 
positive. A similar situation occurs in the present case. 
For the state ([3]), we find 



( a tV) = ^n(n-l)/?. 

n=0 

oo 

(o+a) = np„ ; » 



3|g| 4 (3 + 2|£| 2 ) 

" ^(i - m r/2 1 

ie| 2 (l + 2|e| 2 ) 



iV2(l-|C| 2 ) 5 / 2 ' 



(17) 



where the normalization constant is given by 



N- 



(i-|£l 2 ) 3/2 



(18) 



From (fTT)) . we find that Q becomes negative for |£| < 
0.43, which is satisfied for r < 0.46. We emphasize that 
the Wigner function has negative region for all values of r, 
and thus the photon-subtracted squeezed state is nonclas- 
sical for all r, though it does not exhibit sub-Poissonian 
photon statistics above certain squeezing threshold. 



III. MODELS OF DECOHERENCE 

We next consider how this state evolves under deco- 
herence. The decoherence of the single-mode state §S§ 
can be due to decay of photons to the reservoir or due to 
phase damping. 

A. Amplitude decay model 

When the photons decay to reservoir, the correspond- 
ing Markovian dynamics of the state is well described by 
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the following equation: 

—p = —nicfiap — 2apa) + pa) a) , (19) 

where k is the rate of decay. The solution of this equation 
can be written as 



p(t) = Y,Pn,n'(t)\n){n'\ 



(20) 



where the density matrix element p n>n i (t) can be found 
by using the Laplace transformation and the iteration 
methods [13]. To see this, let us start with the time- 
dependent equation for p n ,n' '■ 



Pn,n> = ~ K { n + n')pn t n> + 2ny/ (n + l)(n' + l)p n +l,n' + l ■ 

( 21 ) 

Using the new subscripts q = n — n' and p = (n + n')/2, 
the above equation transforms into 



Pp, 9 = -2K P p p<g + 2n v / (p+l) 2 - (q/2) 2 p p+hq . (22) 

Taking Laplace transformation of (|22|) and using the orig- 
inal subscript n and n', we can write the time-dependent 
solution for the density matrix elements as 



Pn 



,(t) = e -«*(«+»') ^ J( n+r C r ) ( n '+ r C r ) 



x(l-e- 2Kt ) r p n+r , n , +r (t = 0) , (23) 

where for the single-photon subtracted squeezed vacuum 

1 (£/2)<™+ r+1 )/ 2 (£* /2)(™'+ r+1 )/ 2 
Pn+r,n'+r(t = 0) = - n+r+1 ^ ^ »/ +r+1 y 

(n + r + l)!(n' +r + 1)! 

x- ■ (24) 

y/{n + r)\(n' + r)\ 

We next calculate the parameter Q [Eq. (fT6|) ] for the 
state (|2T)j) . We have found that 



(a t2 a 2 ) = y^n(n - l)p n , w (t) , 

ra=0 
oo 

(a* a) = np n ^ n (t) , 



(25) 



n=0 



where p n ,n(t) is given by (I23|) and (f24|) for n = n' in case 
of state \ip)i. Using Eqs. (j25"|) . we plot Q with time in 
Fig. O It is easy to see that at long times (nt — » oo), Q 
vanishes. This is because at this limit, p n ,n{t) vanishes 
for all non-zero n and po.o(t — ► oo) = 1, i.e., the state 
decays to vacuum. Thus the averages (|25*|) vanish and Q 
also vanishes. 



_Z. Evolution of Wigner function 

The evolution of the Wigner function is governed by 
the following equation: 

d 2 



dW 



d 



d 



da a da* a 



dada* 



W(a,a*) . (26) 



The solution can be written as 
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W(a,a*,t) = 



tt(1 - e- 2fit ) 

\a ~ a a e~ Kt \ 2 
(1 - e- 2 «*) 



d 2 Q lf(a„,aJ,0) 



cxp < —2 



(27) 



where W^cko, (Jq, 0) is the Wigner function of the initial 
state. It is easy to verify this solution putting Q27p in 
(|2"6"| . The time-evolution of the Wigner function of the 
squeezed vacuum state \tp) can be easily calculated ana- 
lytically using the following integral identity [HI, [l^| : 



^Jt 2 - Apv 



exp 



p 2 1/ o z a ZQ 

—a a — ■ 

T T 

pz 2 + vz* 2 + t\z\ 2 



d aexpf— |a| ] exp a a* =■ H ■= 

T T JT Jt 



Apv 



(28) 



Using Eq. (fT5|) and the above identity in Eq. (|2"7|) . we 
get the following: 



W(a,a*,t) 







2|q| 2 


4 


exp 


l_ e -ait« 



7r(l-e- 2Kt ) y'a 2 - 4|c| 
"6 2 c* +b* 2 c + a\b\ 2 



x exp 



where 



a 2 - 4 c 



-2k* 



(29) 



2 cosh(r) + 2 
2a*e~ Kt 



1 - e- 2Kt 

i0 



1 - e- 



c = e 



sinh(r) 



(30) 



Clearly the Wigner function of the squeezed vacuum state 
is Gaussian at all times. 

We now calculate the time-dependence of the Wigner 
function of the state The initial Wigner function, 

as given by ([14]) . can be rewritten as 



W(a,a*,0) = -D 

7T 



\=2 d\ 



(31) 

Using (|31[) and (j2"T|) , we can find the following expression 
for the Wigner function: 

W(a,a*,t) = f~) l _l_ 2Kt D J rf 2 a exp[-A|a | 2 ] 



x exp 



a — age 



1-e 



-2/ct 



(32) 



A=2 



where <So = ao cosh(r) — a$e %e sinh(r). Simplifying the 
above expression using Eq. (|28|) . we get 



W(a,a*,t) = 



32P 



exp 



2|q| 2 
T-e-^* 



7r(l-e- 2 «*) (a 2 - 4|c| 2 ) 5 / 2 



x exp 



6 2 c* + b* 2 c + a\b\ 2 



a 2 - 4 c 



, (33) 
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where 



(1 -£ 2 ){sinh(2r)(e i( V- 
+2[(x + l) 2 +4xsinh 2 (r)]} 
+2{(1 + x 2 ) cosh(2r) + 2x}\b{< 



181*2 



b* 2 ) 



and 



1 



(34) 



(35) 



Clearly, the Wigner function is non-Gaussian due to the 
presence of the polynomial P. This becomes negative 
when the polynomial P becomes negative. In Fig. [6j 
we have plotted C'(a, a*,t) = P/(a 2 - 4|c| 2 ) 5/2 in phase 
space to show the negative region for Wigner function. 

Note that at the center of the phase space (a = a* = 
0), the Wigner function is maximally negative. At the 
center, 



C'(0,0,t) 



2(1 - x 2 ){{x + l) 2 + 4xsinh 2 (r)} 
(a 2 -4|c| 2 ) 5 / 2 



(36) 



which becomes negative when (1 — x 2 ) becomes negative. 
This leads to the following condition: 



nt < nt 



1 







ln(2) 



(37) 



which is independent of the squeezing parameter r. Thus 
the Wigner function has certain negative region for the 
time t < to = ln(2)/2ft. However the negative value does 
depend on the squeezing parameter. 

We hnd from Eq. (fM)) that at t = to (i.e, when x 2 = 1), 
P becomes a circle in phase space. Thus beyond t > to, 
the ellipse 16C"/(1 — e~ 2Kt ) = constant interchanges its 
minor and major axes. We show this behavior in Figs. [7| 
for different values of r. Note that at times much larger 
than decoherence time-scale l/n (i.e., for nt — > 00), P — > 
2 and thus becomes constant throughout the phase space. 

Using Eq. (|3"3"|) , we show the variation of Wigner func- 
tion at different time-scales in Figs. HJ It is easy to see 
how the negative region of the Wigner function gradually 
diminishes. At long times nt — > 00, the Wigner function 
becomes 



W{a, a*, 00) = — e 

7T 



-2|a| 2 



(38) 



which corresponds to vacuum state. We have shown this 
in Fig. [IJb). This can also be understood from Eq. (f!?3")) . 
For nt — * 00, po,o approaches unity, whereas all other 
density matrix elements vanish. This means that at long 
times, the state decays to vacuum, as we have discussed 
earlier. 

We next study the time-evolution of the Wigner func- 
tion for the case of large squeezing, i.e., large values of 
Q. In this case the single photon subtracted squeezed 
state becomes similar to a Schrodinger cat state. For 
large times, such an optical cat state decays to vacuum. 
Thus the Wigner function becomes Gaussian, as dis- 
cussed above. We show this evolution for large squeezing 
in Figs. M 



B. Effect of phase damping 

We now study the effect of phase-damping on the state 
Such damping can be described by the following 
master equation: 



P 



-K p (A^Ap - 2A P A^ + pA^A) , 



(39) 



where A = a^a is the number operator and k p is the 
corresponding rate of decoherence. The solution of this 
equation can be easily found as (|20[) where 



Pn,n'(t) = 6Xp[-(rt - ri) 2 K p t\p ntn ,(0) 



(40) 



It is easy to see that only the diagonal elements p nj „ do 
not decay due to dephasing. Thus at long times, we can 
write 



Pit 



00 



ra=0 



.(0)|n>M 



(41) 



which refers to a mixed state. 

Using Eqs. P5|l and ([4"0|) , we next calculate the pa- 
rameter Q. We find that the averages (a^a 2 ) and (a^a) 
do not depend upon time, because in case of phase damp- 
ing p nt n(t) = Pn,n(0). Thus Q remains the same for all 
times. However, the corresponding Wigner function has 
certain time-dependence. We find that at long times, the 
Wigner function becomes 



W(a, a* ,00) 



n=0 



(42) 



where W\ n \/ n \(a,a*) is the Wigner function of a Fock 
state \n) as given by 



W ln){nl (a,a*) = (-l)"-e- 2 HX(4|a| 2 ) 



(43) 



The function (|4^|) refers to a highly nonclassical state. 
It is interesting to note that all the Fock states have 
independent contributions to the Wigner function at long 
times, weighted by their initial population p n>n (§). On 
the other hand, in case of decoherence due to photon- 
number decay, only the vacuum state survives. In Fig. 
[TOl we plot the Wigner function (|42f in phase space for 
different squeezing. Note that the Wigner function has 
negative region at long times representing nonclassicality 
for all r, even if the state does not exhibit sub-Poissonian 
statistics for r > 0.46 (because Q is positive). In fact, if Q 
is positive, it does not mean that the the state is classical. 
In such cases, we have to use other parameters to test the 
nonclassicality. Several parameters have been introduced 
in this context [H, [2(| ■ We can use hierarchy of these 
parameters which have been shown to be especially useful 
in context of cat states. Here we illustrate the utility of 
one such parameter, e.g., the A3 parameter as defined by 

M 



A, = 



detfm (3) l 



det[^( 3 )] -det[m( 3 )] ' 



(44) 
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FIG. 10: Wigner function in phase space at long times in presence of phase damping for (a) r = 0.31 and (b) r = 0.8. 



where 

(1 mi m 2 \ I 1 Mi M2 \ 

mi m 2 m 3 J , m (3) = Mi M2 M3 I , 
m 2 m 3 m 4 y y ti 2 M3 / 

(45) 

m s = (a^ s a s ), fi s = ((ata) s ) and det indicates determi- 
nant of the matrix. The state exhibits phase-insensitive 
nonclassical properties if A3 lies between and -1 
For the state we have found that A3 remains nega- 
tive for |£| < 0.6 which corresponds to r < 0.7. Clearly 
A3 is a stronger measure of nonclassicality than Q be- 
cause it leads to a larger upper bound of r to exhibit 
nonclassicality. Further, comparing the Wigner functions 
in Figs. [10] with those at t = [see Figs. [2], we find that 
the Wigner function varies very slowly with time for small 
squeezing. But for large squeezing, the variation is faster. 
Although we can conclude that phase damping leads to 
much slower decoherence than amplitude damping. 

IV. CONCLUSIONS 

In conclusion, we have studied how a class of non- 
Gaussian states evolves in presence of decoherence. We 
have considered a single-photon subtracted squeezed vac- 
uum state, the Wigner function of which is similar to 



that of a Schrodinger kitten state. We have found an up- 
per bound for squeezing parameter for which this state 
exhibits sub-Poissonian photon statistics. However, the 
state remains nonclassical for all values of the squeezing 
parameter because the Wigner function becomes nega- 
tive around central region in phase space. Next, we have 
studied how the state evolves in presence of two different 
kinds of decoherence, viz., amplitude decay and phase 
damping. We have found analytical expressions for the 
time-evolution of the state and the Wigner function in 
both cases. In case of amplitude decay, the Wigner func- 
tion loses its non-Gaussian nature and becomes Gaussian 
at long times, corresponding to vacuum. On the other 
hand, phase damping leads to much slower decoherence 
than amplitude damping. The state remains nonclassical 
at long times. 
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